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Abstract 

The effects of quantum lattice fluctuations on the Peierls transition are stud- 
ied within the one-dimensional Holstein molecular crystal model by means of 
exact diagonalization methods. Applying a very efficient variational Lanczos 
technique, the ground-state phase diagram is obtained in excellent agreement 
with predictions of recent density matrix renormalization group calculations. 
The transition to the charge-density-wave regime is signaled by a strong 
increase in the charge structure factor. In the metallic regime, the non- 
universal Luttinger liquid parameters (charge velocity and coupling constant) 
are deduced from a finite-size scaling analysis. The variational results are 
supported by a complete numerical solution of the quantum phonon Holstein 
model on small clusters, which is based on a well-controlled phonon Hilbert 
space truncation procedure. The metallic and charge-density-wave phases 
are characterized by significant differences in the calculated optical absorp- 
tion spectra. 

PACS number(s): 71.38.+i, 71.45.Lr, 72.10.Di 
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I. INTRODUCTION 



Many quasi one-dimensional (ID) materials, such as the organic conjugated polymers 
[e.g., (CH) Z ] and charge transfer salts [e.g., TTF(TCNQ)] or the inorganic blue bronzes [e.g., 
K0.3M0O3] and MX-chains |T[, undergo a Peierls instability in the half-filled band case at 
temperatures between 50 K and 250 K, driven by the electron-phonon (EP) interaction 0. 
Most theoretical studies of these systems concentrate on the ID SSH [[J and Holstein ||] 
models, where the phonons interact with the electrons by modifying the electron hopping 
matrix element and on-site potential, respectively. Frequently the lattice degrees of free- 
dom were treated classically. However, it has been argued that for most quasi-lD systems 
the lattice zero-point motion is comparable to the Peierls lattice distortion, which makes 
the rigid lattice approximation questionable || . Although the problem, whether the dimer- 
ized ground state survives the quantum phonon fluctuations, has been addressed by several 
numerical |||| and analytical P-pH approaches, it has to date resisted a complete theoret- 
ical solution. Moreover, lattice dynamical (quantum phonon) effects should be included in 
any theoretical analysis of the extraordinary transport and optical phenomena observed in 
Peierls-distorted systems [jT2|JT~3[1 . The discussion of optical properties, e.g., of the optical 



absorption, poses, however, an extremely complicated many-body problem that cannot be 
solved analytically without further approximations []13} , |n)| . Perhaps, at present, the most 
reliable results for models of electrons strongly interacting with quantum phonons come from 
finite-cluster calculations supplemented by a careful finite-size analysis. 

Encouraged by this situation, in this work we present a purely numerical (exact diagonal- 
ization) study of the ID spinless fermion Holstein model ||16|| , with a view to understanding 
the effect of quantum lattice fluctuations on the Peierls dimerization and optical absorption 
spectra in both the metallic [Luttinger liquid (LL)] and insulating [charge-density-wave 
(CDW)] phases. In a particle-hole symmetric notation the Holstein Hamiltonian reads 



H = -t£(4c m + cj +1 cj - yV^o^&l + b t ) (n t - 1) + Tiu^iAk + ±) , (1) 



where c[ (b\ ) are the electron (phonon) annihilation [creation] operators, and n i = cfcj. In 
Eq. (|IJ), the free-electron transfer amplitude t is restricted to nearest-neighbour hopping, a 
dispersionsless Einstein phonon uj(q) = ujq is coupled to the local electron density, and the 
phonons are treated within harmonic approximation. In the atomic limit (t = 0), e p gives 
the well-known Lang-Firsov polaron binding energy. Rescaling 7i — > 7i/t and measuring 
all energies in units of t, it is convenient to introduce the adiabaticity parameter a = hu^/t 



and two dimensionless EP coupling constants, A = e p /2t and g = ^Je p /hujQ, in order to 
characterize the weak (A <C 1) and strong coupling (A ^> 1 and g 3> 1) situations in the 
adiabatic (a <C 1) and anti-adiabatic (a ^> 1) regimes. For the single-carrier case, the 
Holstein Hamiltonian was extensively studied in the context of the polaron problem. A 
practically complete numerical solution, including ground state and spectral properties, is 
now available from ED (exact diagonalization) and DMRG (density matrix renormalization 
group) calculations JT7Hl9|. At half-filling, a number of different analytical and numerical 
methods, including strong coupling expansions J7j , variational approaches ]20| and renormal- 
ization group arguments plfl , as well as world-line quantum Monte Carlo (WL QMC) [0, 
Green's function Monte Carlo (GF MC) || and DMRG techniques [p^j , have been used to 
determine the phase boundary between metallic and insulating behaviour. 

In this paper, we follow our strategy pursued in recent work on the self-trapping prob- 
lem in the single-electron Holstein model |[23| , and apply two different numerical meth- 
ods: (i) a complete exact diagonalization (ED) of the Holstein model preserving the full 
dynamics and quantum nature of phonons and (ii) a variational Lanczos scheme based 
on the inhomogeneous modified variational Lang-Firsov transformation (IMVLF) ||24|| . It 
is natural for the first method to be limited to rather small clusters (with sites) and 
moderate EP coupling strengths (A, g), mainly due to the necessity of a phonon Hilbert 
space truncation (retaining at most M phonons; for details of the numerical method see 



Ref. [25]). By combining ED with Chebyshev recursion and maximum entropy meth- 



ods [2c,2q|, we are able to discuss the dynamical properties of the spinless fermion Hol- 
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stein model, such as the optical conductivity (see Sec. IV). On the other hand, using the 
second method, we can study the ground-state properties of fairly large systems (Sec. II), 
which enables us to carry out a finite-size scaling (cf. Sec. III). Within the IMVLF-Lanczos 
approach, we treat the phonon subsystem by performing first of all a canonical transfor- 
mation, H = WHU, U = e^^'le^We- 5 ^ 1 ), where <Si(A) = -^EiA(b\ - k), 
£2(7,7) = -9Y,i{b\ - &<) (7 + 7«i), and 5 3 (r) = |lnrEi(&J&l - W are designed to de- 
scribe static displacement field (Ai), non-adiabatic polaron (7,7), and squeezing (r) effects, 
respectively. Next, we approximate the eigenstates of the transformed Hamiltonian by 
the variational product states |$v) = \& p h) ® |^ei) an d average H, over the phonon vacuum, 
H = (^phl^l^pft), which leads to an effective electronic Hamiltonian 

H = f a tf + 7)^ + /a[ 7 2 - 7 + 27(7 - 1)] - e_9272r2 EfoVi + 

i i 

-(l-T)EAK-l) +t/EA + ^(- 2 + r- 2 ) + T^EA 2 - (2) 

Here 77 is a Lagrange multiplier ensuring the constraint Yli A i — 0. Employing the Hellmann- 
Feynman theorem, the iV + 3 variational parameters are obtained by iteratively solving the 
extremal equations for the corresponding energy functional Eo({Ai}, 7, 7, t 2 ) in combination 
with the Lanczos recursion algorithm. 

II. PHASE DIAGRAM 

Previous results for the ground-state phase diagram of the Holstein model at half-filling 
obtained by WL QMC and GF MC || simulations showed significant discrepancies in 
the region of small a (0 < a ^ 1). Only very recently Bursill et al. provided more 
reliable information from level crossings in their DMRG data. Applying, in a first step, our 
variational Lanczos scheme, we consider the effective model (0) on chains of even length 
with up to 16 sites and periodic ( ant iper iodic) boundary conditions if there is an odd (even) 
number of fermions in the system. To elude the problem of trapping in metastable minima 
of the energy functional Eq, we start the variational Lanczos iteration with different initial 
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configurations {Aj, 7, 7, r 2 }, close to the metallic or the dimerized phase. As one might 
expect, in the dimerized phase the iteration always converges to a ground-state with stag- 
gered dimerization Aj = A(— l) 4 , while in the metallic phase Aj = 0. Thus, at half-filling, 
our inhomogeneous variational wave function becomes exactly the staggered (SMVLF) one 
used in Ref. [2~0"|, and the effective Hamiltonian 7Y(A,7,r 2 ) can be solved easily, also for 



the infinite system. However, the phase diagram given in |20j for the infinite system is only 
tentative, and for small a the determination of g c is not clear. Moreover the infinite system 
is never really gapless within the variational approach, because A remains nonzero, although 
it becomes very small for weak EP coupling. This situation changes as soon as the system 
is finite. The dimerization A then switches from zero to a finite value at a critical coupling 
g c (a,N), where g c is nearly independend from the system-size A" for large a (a ^ 1), but 
decreases with A" for small a. 

Proceeding this way we get the IMVLF transition lines g c {ct, N) depicted in Fig. 1. Most 
notably we found that the IMVLF phase boundary, separating metallic (LL) and insulating 
(CDW) phases, in the whole parameter regime, agrees surprisingly well with the very recent 
DMRG results []2I| (open squares). Also in the phase diagram are the transition points 
obtained by WL QMC § and GF MC |. 

In a second step, performing exact diagonalizations of the full Hamiltonian for systems 
with 6 and 10 sites and up to 30 phonons, we want to demonstrate that the phase boundary 
determined from the effective model (0) is consistent with what is obtained for the quantum 
phonon Holstein model (|]). To this end, we have calculated the static charge structure 
factor, 

XW = 7^E^ (R ^ ) (^), (3) 

shown in Fig. 2 as a function of g at low (a), intermediate (b), and high (c) phonon fre- 
quencies for both, the variational and the exact solution. Increasing the EP coupling at 
fixed phonon frequency, the smooth variation of xi^) m t ne metallic phase is followed by a 
strong enhancement at about g^ MWLF \ unambiguously indicating the formation of a CDW. 



The discontinuous jump-like behaviour of ;\;( IMVLF ) (7i~) at g c is an apparent shortcoming of 
the variational approach. This resolves also the open question in [2(J, whether the two- 
minimum structure of the variational solution at large a is an artifact. For a — > and 
N — > oo, where the IMVLF approach becomes exact, we found a continuous crossover in 
X (IMVLF) (vr) as well. 

Obviously, in the CDW-like phase, a larger number of phonons (M) is required to achieve 
a satisfactory convergence of the ED data (see Fig. 2 (b)). Furthermore, it is interesting 
to compare the behaviour of the kinetic energy (Ek in ), given by the average of the first 
term of ([[]), in the adiabatic, non-adiabatic and anti-adiabatic regimes (see insets). For 
the adiabatic case, the kinetic energy is only weakly reduced from its noninteracting value 
((Ekin) = —4) in the metallic phase. By contrast, in the anti-adiabatic regime, we observe a 
strong reduction of (Ekm), which can be attributed to the formation of a strongly correlated 
polaronic metal below the CDW transition point. 

Coming back to the phase diagram, in the adiabatic regime, our results seem to confirm 
that there is no long-range order for sufficiently small EP coupling, which is consistent with 



the predictions of Refs. p7|,[8|. At a = 0, the critical coupling converges to zero, as expected 
for the adiabatic Hamiltonian (M — > oo; 7 = 0, r 2 = 1). In the regime < a ~ 1, however, 
the precise determination of g c is somewhat difficult. Maybe the discrepancy between the 
predictions of IMVLF, GF MC, DMRG on one side and WL MC on the other side, of how 
the critical A c scales to zero with a — > 0, results from this ambiguity. 

In the strong-coupling non-adiabatic regime (g 2 , A ^> 1), the results of the different 
numerical approaches approximately agree. Also the analytical approach, giving the exactly 
soluble XXZ model [fF|,p8| within second order perturbation theory (with respect to t) 

H xxz = j(2a - g 2 a - V 2 ) - e"^ £ ((S+Sr +1 + SrSf +1 ) - V 2 e° 2 S*S* +1 ) (4) 

i 

with 

V n {a,g) = 2^ — , (5) 

a s _^ ss. 



works very well. The (Kosterlitz-Thouless) phase transition line is given by the condition 
V 2 {a, g 2 )e g2 /2 = 1 (dashed curve in Fig. 1). For a — > oo (anti-adiabatic limit) there is no 
dimerization if A is finite. 

To get a feeling about the accuracy of the different analytical and numerical techniques, 
we have compared in Fig. 3 the ground-state energies at high phonon frequencies, where the 
small polaron approximation is justified. For the adiabaticity ratio a = 10, used in Fig. 3, 
the transition to the CDW phase takes place at about g c ~ 2. 

Restricting ourselves to the metallic phase, the ground-state energy of the XXZ model 
is easy to evaluate. It is (per site) 

with cos[/i] = V^e 9 /2. Alternatively, calculating the polaron self-energy of the Holstein 
model (0) within standard second-order (Rayleigh-Schrodinger) strong-coupling perturba- 
tion theory (SCPT) and omitting the residual polaron interaction, the small polaron band 
dispersion becomes JTS 

2 

^SCPT = _9_^_ + y 2 _ 2e -g- cog R _ Vi e -g* ^ 2R _ (?) 

Then, for the half-filled band case, the ground-state energy (per site) takes the form 

£o SCPT = + | , (8) 
where the (first-order) polaron ground-state energy, 

can be obtained from (^]) by setting A± — 0, 7 = 1, 7 = —1/2, and r 2 = 1. 

Both E^ xz and Eq CPT are also depicted in Fig. 3, and in order to visualize the higher- 
order corrections we have shifted all energies by the standard small polaron term ^q 1 '*. First 
of all we see that the higher-order corrections, originated by the residual polaron-phonon 
interaction, are most important at intermediate couplings g ~ 1, where the polaron band 



structure significantly deviates from a rescaled cosine tight-binding band [fT8 1 . As already 
mentioned above, the IMVLF-Lanczos results, extrapolated to iV = oo, coincide with the 
variational SMVLF solution |20[ , but give higher ground-state energies than the XXZ model 
in the intermediate-to-strong coupling regime. This is because the non-adiabatic polaron 
effects are only included to the lowest order of approximation (remind that the model (^) 
was obtained performing the average over the zero-phonon state). Including second order 
corrections by SCPT, we found a much better agreement with the exact data (full circles). 
Notice the substantial lowering of the energy with respect to the first order result at large 
EP couplings [g — 2 implies A = 20 (!)], which results from the momentum independent 
shift in the kinetic energy [jny^J. As expected, the finite-size effects due to the lattice 
discreteness are most pronounced in the weak-coupling regime. 



III. LUTTINGER LIQUID PARAMETERS 



According to Haldane's Luttinger liquid conjecture |30|], ID gapless systems of interacting 



fermions should belong to the same universality class as the Tomonaga-Luttinger model. As 
stated above, the Holstein system is gapless for small enough coupling g. Thus it is obvious 
to prove, following the lines of approach to the problem by McKenzie et al. |y, whether our 
IMVLF-Lanczos data shows a finite-size scaling like a Luttinger liquid. 

For a LL of spinless fermions, the ground-state energy Eq(N) of a finite system of N 



sites scales to leading order as [31 



E (N) = tiu p 
N e °° 6A 2 



(10) 



where denotes the ground-state energy per site for the infinite system and u p is the 
velocity of the charge excitations. If E±i(N) is the ground-state energy with ±1 fermions 
away from half filling, to leading order the scaling should be 

7T?/ 

E ±1 (N) - E (N) = ^ . (11) 
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K p is the renormalized effective coupling (stiffness) constant (for a more detailed discussion 
see Ref. [g). 

As becomes evident from Fig. 4, our IMVLF-Lanczos data match both scaling relations 
for all regimes of parameters a and g with great accuracy. Let us stress that the IMVLF- 
energies for E±i(N) correspond to a inhomogeneous solution for the displacement fields 
Ai, which deviates from a uniform or staggered ordering, and therefore cannot be obtained 



within the simple SMVLF scheme |20 



In the plot of the LL parameters, Fig. 5, increasing error bars at the end of the curves 
indicate the phase transition to the CDW state. Here the scaling relations ( |10D and fllT]) 
no longer hold. Perhaps surprising, the velocity of charge excitations u p , reflecting the 
behaviour of the kinetic energy (cf. Fig. 2 insets), agrees fairly well with the results of 
McKenzie et al. ||, while the coupling constant K p is always bigger than one, indicating an 
attractive interaction. We believe this is an artifact of our variational treatment. For the 



XXZ model K p — ► 1/2, i.e. the phase transition is of infinite order with a Kosterlitz- 
Thouless order parameter ~ e~ l ^ 9 ~ 9c \ Note that the parameters u p and K p we show, are 
of course those for the effective Hamiltonian (|2|) , which coincide with the parameters of the 
true Holstein model only more or less. 

Unfortunately, extracting a similar scaling behaviour from our ED data seems to be 
extremely complicated, mainly because two finite-size dependences, those with respect to 
the system size N and maximum phonon number M, are mixed. Moreover, the memory 
limitations of the present day parallel computers impose severe restrictions on the lattice 
sizes (N < 10), that can be treated by ED with adequate accuracy. 



IV. OPTICAL RESPONSE 

One of the physical quantities which contains extremely valuable information about the 
low-energy excitations in polaronic metals and CDW systems is the optical conductivity, 
c(uj), usually determined from reflectivity measurements. As will be shown in the following, 
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the optical absorption spectra of nearly free electrons, small polarons and CDW insulators 
differ essentially. 

The real part of o~(u) contains two contributions, the familiar (coherent) Drude part at 
uj = and a so-called "regular term", a re9 (uj), due to finite-frequency dissipative optical 
transitions to excited quasiparticle states. In spectral representation (T = 0), the regular 
part takes the form |25|j3^1 



= E IWE ' (c k +1 C ^»^» 2 8[u - (E m - E )} , (12) 

m>0 ~ ^0 



where a re9 {u) is given in units of 7re 2 and we have omitted an 1/N prefactor. In (|12|). 



the summation is taken over the complete set of eigenstates \& m ) with excitation energies 
uj = (E m — E ) in the subspace of N/2 (spinless) fermions (half-filling). For the discussion of 
the optical properties it is useful to consider also the u-integrated spectral weight function 

S re9 (u) = I" duj'a re9 (J). (13) 
Jo 



The evaluation of dynamical correlation functions, such as (|12|), can be carried out by 
means of very efficient and numerically stable Chebyshev recursion and maximum entropy 
algorithms Nevertheless, due to the huge size of the Hilbert space in this EP model, 

we are currently restricted to a lattice size of 6 sites (M = 30; periodic boundary conditions) 
if we want to calculate the conductivity in a wide range of EP coupling strengths. 

Typical optical absorption spectra for the ID half-filled Holstein model of spinless 
fermions are given by Figs. 6, 7, and 8, in the metallic (a) and CDW (b) phases, at charac- 
teristic phonon frequencies corresponding to the adiabatic, intermediate, and anti-adiabatic 
regimes, respectively. 

For low phonon frequencies (a < 1) and weak EP couplings (see Fig. 6 a), the peak 
structure may be easily understood in connection with the non-interacting tight-binding 
band dispersion E$ = —2 003^, where the allowed K values are K — 0, ±7r/3, ±27r/3, 
and 7r for a six site system. Obviously, we found the first transitions with non-negligible 
(electronic) spectral weight (cf. S reg (uj)) at frequencies that approximately correspond to 
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to the discrete free electron Bloch states of our finite system and its vibrational satellites. 
Accordingly the first and second group of excitations originate from transitions where the 
momentum of one electron is changed from ±7r/3 to ±27r/3 and ±7r/3 to ±7r, respectively. 
Note that in ([12]) an optical transition only takes place within the K — sector (l^o) 
carries K = for the half-filled band case). Thus a phonon with opposite momentum must 
be absorbed in order to ensure momentum conservation during a single-particle excitation 
process. Of course, in the Holstein model, K is the total momentum of the coupled EP 
system, and, at any g > 0, there is a finite overlap of the ground state with all the excited 
states belonging to the same K sector. The most relevant point is that in the metallic 
phase the absorption threshold should tend to zero as the number of sites increases, i.e., the 
low-energy (finite-size) gap vanishes. 

The optical absorption spectrum in the strong EP coupling regime is quite different from 
that in the LL phase. It can be interpreted in terms of strong electron-phonon correlations 
and corroborates the CDW picture. For g > g c the electronic band structure is gapped (at 
the edge points K = ±7r/2), and we expect that now the low-energy gap feature, observed 
in Fig. 6 (b), survives in the thermodynamic limit N — > oo. Unfortunately finite-size effects 
prevent a precise extraction of the CDW gap from our optical ED data. The broad optical 
absorption band found above this gap is produced by a single-particle excitation accompanied 
by multi-phonon absorptions and is basically related to the lowest unoccupied state of the 
upper band of the CDW insulator. The lineshape reflects the phonon distribution in the 
ground state (see Fig. 9). The most striking feature is the strong increase of the spectral 
weight contained in the incoherent part of optical conductivity. This becomes evident by 
comparing the magnitude of S re9 (oo) in the weak and strong coupling situations. Moreover, 
employing the f-sum rule for the optical conductivity and taking into account the 
behaviour of the kinetic energy as function of g (see Fig. 3), we found that in the metallic 
and CDW phases nearly all the spectral weight is contained in the coherent (Drude) and 
incoherent (regular) part of Re o~(uj), respectively. That is, in the CDW state the transport 
is dominated by inelastic scattering processes. 
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In the non-adiabatic region a ~ 1, where the phonon frequency becomes comparable to 
the electronic bandwidth (level spacing), the situation is not much different (see Fig. 7). 
Again, in the weakly interacting case, the optical absorption can be understood in terms of 
electronic transitions within a tight-binding band and phonon satellites. Crossing the CDW 
transition point, a pronounced redistribution of spectral weight from the Drude to the regular 
part of <t{uj) is observed. As can be seen by comparing Figs. 6-8 and Fig. 9, the heights 
of the jumps in the ^-integrated conductivity, being directly related to the probability of 
the corresponding m-phonon absorption processes, give a measure of the weights of the 
m-phonon states in the ground state. 

Finally, we consider the optical response in the anti-adiabatic regime. For weak inter- 
action the picture remains the same as in the above adiabatic and intermediate cases (see 
Fig. 8). If we increase the EP coupling g, the electrons will be heavily dressed by the phonons 
(which now can follow the electron instantaneously) and the formation of less mobile small 
polarons takes place (cf. (E kin (g)) and u p (g) shown in Fig. 3 (c) and Fig. 5 (a), respectively). 
As a consequence the coherent transport becomes strongly suppressed in the LL phase. Since 
the (renormalized) coherent bandwidth of the polaron band is rather small, the finite-size 
gaps in the band structure are reduced as well, and the CDW gap (A C dw ~ 2A) may be 
identified with the optical absorption threshold. Again we found that in the CDW phase 
multi-phonon absorption processes dominate the optical response. 

V. SUMMARY 

In this paper we have studied the Peierls instability and the optical absorption in the 
half-filled spinless fermion Holstein model by means of finite-lattice diagonalizations. We 
have shown that the simple variational (IMVLF) Lanczos approach can be successfully used 
to determine the ground-state properties of the Holstein model, in particular the phase 
diagram. The calculation of the optical properties did require a complete diagonalization of 
the model, preserving the full dynamics of the phonons. 
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Our results confirm previous findings that at sufficiently weak electron-phonon (EP) 
interaction the system resides in a metallic (gapless) phase, described by two Luttinger 
liquid parameters. The renormalized charge velocity (u p ) and correlation exponent (K p ) 
were obtained from finite-size scaling relations, fulfilled with great accuracy Increasing 
the EP coupling, the system undergoes a Peierls transition to an insulating (gaped) phase, 
reflected in a strong increase of the charge-density correlations with momentum n. The 
crossover between Luttinger liquid and charge-density-density (CDW) behaviour is found 
in good agreement with exact diagonalization and density matrix renormalization results for 
the quantum phonon model. In the non-adiabatic strong-coupling limit, where the charge 
carriers are polaronic, the IMVLF-Lanczos phase boundary lies close to the analytic findings 
for the XXZ (small polaron) model. The transition to the CDW state is accompanied by 
significant changes in the optical response of the system. Most notably seems to be the 
substantial spectral weight transfer from the Drude to the regular (incoherent) part of the 
optical conductivity, indicating the increasing importance of inelastic scattering processes 
in the CDW (Peierls distorted) regime. 

The numerical results made clear that a dynamical treatment of the lattice degrees is 
necessary in the intermediate frequency and coupling region because the energy scales are not 
well separated. This should be a matter of relative importance modeling the inorganic spin 
Peierls materials |[34| , e.g. CuGe03, where the spin exchange interaction and the relevant 
phonon frequencies are of the same order. Indeed, applying our numerical techniques to a 
frustrated Heisenberg spin— ~ chain with dynamic spin-phonon coupling shows that in the 
non-adiabatic regime the spin-Peierls transition takes place at about g c ~ 1, which implies 
for CuGeOa a intermediate to strong coupling situation . 
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FIGURES 



FIG. 1. Ground-state phase diagram of the ID Holstein model of spinless fermions at 
half filling, showing the boundary between the Luttinger liquid (LL) and charge-density- 
wave (CDW) states. The IMVLF-Lanczos results are compared with the predictions of 
different analytical and numerical approaches. 

FIG. 2. Charge structure factor x(n) and kinetic energy (E kin ) (inset) as a function of 
the EP coupling g in the adiabatic (a), non-adiabatic (b), and anti-adiabatic (c) regimes. 

FIG. 3. Dependence on the EP coupling of the ground state energy (shifted by E^). 
The predictions of strong coupling expansions are confronted with IMVLF-Lanczos and ED 
results obtained for a six-site lattice at a — 10. 

FIG. 4. Finite-size scaling of the charge gap E_i(N)—E (N) and the ground state energy 
E (N) (inset) for different values of g (g 2 = 0.1, 0.6, 1.1, 1.6, 2.0 from top to bottom) at 
a = 0.1 (a), 1 (b), and 10 (c). 

FIG. 5. LL parameters u p [charge velocity (a)] and K p [correlation exponent (b)] as a 
function of the EP coupling g. The dashed curves denote the corresponding results for the 
XXZ model at a = 10. 

FIG. 6. Regular part of the optical conductivity a reg (uj) (dotted line) and integrated 
spectral weight S re9 {uj) in the adiabatic weak (a) and strong (b) EP coupling regimes (N = 6, 
M = 30. 

FIG. 7. The same as Fig. 6 but at a = 1. The results in the LL and CDW regimes are 
shown an (a) and (b), respectively. 

FIG. 8. Optical absorption in the anti-adiabatic regime (a = 10) of the ID half-filled 
Holstein model. 
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FIG. 9. Phonon-distribution function \c m \ 2 , as defined in equation (A9) of p5|, shown 
for the parameters used in FIG. 6-8. 
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